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Abstract. The connection between symmetries, time independent and time dependent first 
integrals for autonomous systems of first-order ordinary differential equations is discussed. 
Several new results which have been obtained by applying a reduce package are reported. 


Let x = (xj,... x„) be n functions depending on the variable t and let 

n„(f, x) = x v -a> v (x) = 0 (v=l,...,n) (1) 

be a system of n autonomous differential equations of first order. Throughout it is 
assumed that the o> v are polynomials in x. The infinitesimal generator U of a point 
symmetry is defined by 

U = Z(t,x)d/dt + T la (t,x)d/dx a . (2) 

Summation from 1 to n over twice occurring indices is always understood. The 
infinitesimals £ and r) a have to be determined from the invariance condition U (l) Cl v - 0 
on = 0 where C/ (1) is the first extension of U. The following system of linear partial 
differential equations for £ and the T) a 's is obtained 

dii v /dt+(o a dT) v /dx a - r) a doi v /dx a -a) v (d^/dt +(U a dg/dx a ) = 0. (3) 

The general solution of the system (3) determines the full symmetry group for equations 
(1). This is in most cases impossible. 

Symmetries often make it possible to diminish the number of dependent variables. 
To this end the variables x have to be transformed such that the corresponding 
infinitesimal generator of the symmetry becomes a translation in one of the dependent 
variables. From this it is clear that only those symmetries are useful for the integration 
of (1) for which this transformation may actually be performed. 

Let us now discuss several special solutions of (3) from this point of view: 

(i) For arbitrary co a there exist the obvious solutions 

£ = constant, T 7 „ = 0, t translations (4a) 

£=f(t,x), ri v =f(t,x)cx)„ /arbitrary. (4b) 

These symmetries are useless for the integration of (1) for the reason mentioned above. 
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(ii) If t?„ = 0 for all v and £ does not depend on t, then the system (3) reduces to 

(o a d£(x)/dx a = 0. (5) 

If F(x) is a solution of (5) it is a time independent first integral and the system (1) 
allows the symmetry generator 


U = F(x) d/dt. 


( 6 ) 


Symmetries of this kind are the most useful ones for the integration. 

(iii) Again let i 7 „ = 0 for all v and an exponential time dependence for £ be 
assumed, i.e. 

£(r,x) = e%0c). (7) 

The system (3) reduces to 

*&(*) + &>„ HoM/dx a = 0. (8) 

The solutions of (7) and ( 8 ) are the time dependent first integrals which have been 
proposed by one of the authors (Steeb 1982). If such an integral F(t, x) with exponential 
time dependency can be found, the system ( 1 ) again allows a symmetry generator of 
the form ( 6 ) with a time dependent function F. 

More general symmetry generators may be obtained if no restrictions on £ are 
imposed. This discussion shows the relation between time independent first integrals, 
time dependent first integrals with exponential time dependence and a general symmetry 
generator. 

We now consider several autonomous systems of the form ( 1 ) and determine a 
certain class of their symmetries. It turns out that all calculations in this field are very 
lengthy. For this reason there have been developed two computer algebra packages 
(Schwarz 1984a, b). If polynomial dependency of a predefined degree in x is assumed, 
they determine the time independent and the time dependent first integrals and more 
general infinitesimal generators of a symmetry almost completely automatically. Most 
results described below have been obtained by applying these reduce packages. 

We begin with the Lotka-Volterra equations for three competing populations which 
has been considered by Steeb and Erig (1983) and Steeb et al (1983), namely 

x, =X((1 +ax 2 + bx 3 ), x 2 -x 2 (l-ax,+ / Sx 3 ), x 3 = x 3 (l - bx, - cx 2 ) (9) 

where a, b and c are real parameters. For arbitrary a, b, c there is the first integral 

I(t,x) = {x x +x 2 + x 3 )e0'^ZL^ ( 10 ) 


In addition there is the integral 
_ I(t, x hyJcfxUr e“ (fc+,+m) ' 



( 11 ) 
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if k, l and m are determined from ak-cm= 0, bk + cl = 0. For special values of the 
parameters there exist time dependent integrals of higher order in x. If b= - a, c = a 
there is the third-order integral 

I(t, x) = Of +x\ +x\ +3xi(x 2 +xf) +3x 2 (xf +xf) +3x 3 (xf +x 2 )] e -3 '. (12) 

If b = - a, c = 2a there is the fourth-order integral 

I(t,x) = [x* + xl + x$+ 4x, (x 2 + X 3 ) + 4x 2 (xf + xf) + 4x 3 (xf + xf) 

.2 2 ■ 2 2 . 2 2 \ , . ^ (13) 


+ 6 (xfx 2 +xfxi + X2X3) + 12X[X2X 3 (X2 +x 3 )] e - 
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